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Abstract 

The coupling constant g pal is determined from light cone QCD sum rules. A 
comparison of our result with the ones existing in literature is presented. 
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1 Introduction 



The quark model describes hadrons successfully. In its simplest version mesons are inter- 
preted as pure qq states. Scalar mesons might constitute an exception to this successful 
scheme and indeed their nature is not well established yet [1]. At present there goes great 
deal of investigation about the identity and interpretation of the lowest lying scalar reso- 
nances. One remarkable feature of these particles is that they appear to be rather wide 
[2]-[5]. Other peculiar property of these particles is the observation that they have very 
short lifetimes and large couplings to the hadronic channels such as KK or -kit. It should 
be noted here that the mass of the scalar mesons varies in the region (400 — 1200) MeV 
[6]- ( 

Photoproduction of the p° meson with proton near threshold can successfully be de- 
scribed by a simple one meson exchange, more precisely, by exchange of a scalar a meson. 
An analysis of experimental data requires to know the g pal coupling constant. This coupling 
constant is calculated in frame work of the effective Hamiltonian approach and 3-point sum 
rules in [7, 8], which predict g pal = 2.71 and g pal = 3.2 ±0.6 at m a = 0.5 GeV, respectively. 
Since this coupling constant is one of the important quantities in analyzing the hadronic 
processes involving a meson, its calculation in frame work of the other approaches is needed 
in order to get a reliable conclusion. 

In the present note we study the coupling constant g pay in frame work of the light cone 
QCD sum rules (more about QCD sum rules and its applications can be found in [9, 10]). 
In further analysis we assume that the content of a meson is a pure qq state (although this 
point is still under debate, see for example [11, 12]). 

The coupling of the a meson to the scalar current J°~ = \{uu + dd) can be parametrized 
in terms of a constant A CT : 

(0\J a \a) = \ a . (1) 

In order to determine the coupling constant g pal in frame work of the QCD sum rules, we 
consider the following two point correlator function 

Tln = iJ d 4 xe ip2X (0 |t{ J a {x) J£(0)}| 0) , (2) 

where 7 means external electromagnetic field, and J? and J a are the interpolating currents 
with p and a meson numbers. The physical part of the sum rules can be obtained by 
inserting a complete set of one meson states into the correlator: 

^ (o\.r(x)\a(p 2 )) (p(pi)|^(Q)|o) 

n M = E ZT— ~2 W(P2)\p(Pi)), 2 _ 2 , (3) 
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where p 2 = Pi + q and q is the photon momentum. The matrix element (j}(jp\) Jfc(0) 
Eq. (3) is defined as 
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P) = m p f p e^ , (4) 
where e p is the p meson polarization vector. In general, the (cr\p) matrix can be parametrized 



as 



(*(P2)lA>(Pi)> 7 = e{F 1 ( q 2 )(p iq )e p + F 2 (q 2 )^q)p lp }e^ , 



1 



where e is the photon polarization vector. In this parametrization of the matrix element 
(a\p) , we neglect the terms ~ q^ since q^e^ = 0. From gauge invariance we have 



q^{F 1 ( Pl q)e^ + F 2 (e^q) PllI } = . 
Since photon is real in our case, we immediately get from Eq. (5) 

F 2 (0) = -Fx(0) . 
So, the matrix element (cr\p) takes the following form: 

Hp> 7 = eF(0){( Pl q)e^ - {qe p )p ltl }e^ . 
We can use an alternative parametrization for the pa^ vertex, i.e., 

Cm = — g P *~fd a pP (d a Ap - dpA a ) . 



(5) 



When we compare Eqs. (6) and (7) we see that 



F(0) = 



(6) 



(7) 



(8) 



Using Eqs. (1), (3), (4) and (6), for the phenomenological part of the sum rules we get 



A, 



fp£ v 



P2 -ml pi- m 2 



(PiQ)9iw + Pi 



(9) 



Our next problem is the calculation of correlator function from QCD side. First of all 
we note that the structure of the correlator is such that perturbation contribution to it is 
absent due to the odd number of 7 matrix. So, QCD part of the correlator contains only 
nonperturbative contribution. Calculations lead to the following result for the theoretical 
part of the correlator: 



e(qq) 
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(10) 



where <p and g\ twist-2 and twist-4 photon wave functions, respectively, defined as [13]-[15]: 
(7(9)l^a/3g|0) = ie q (qq) due luqx \ (e a qp - e p q a ) \x<p{u) + x 2 (^i(u) - g 2 ( 



+ 



qx(e a Xf3 - Epx a ) + ex(x a qp - Xpq a ) g 2 {u) \ , 



(11) 



where x is the magnetic susceptibility of the quark condensate, e q is the quark charge. 

The sum rules is obtained by equating the phenomenological and theoretical parts (in 
Eq. (10) it is necessary to perform Fourier transformation first) of the correlator. Perform- 
ing double Borel transformations on the variables p\ = p 2 and p\ = (p + q) 2 on both sides of 
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the correlator in order to suppress the contributions of the higher states and the continuum 
(see [16]— [18]), we get the following sum rules for the g pal coupling constant: 



-(m 2 p /M?+ml/M$) _ ( e u ~ ^d){qq) 

fpX a 



X¥(uo)M 2 E (s /M 2 ) - A 9l (u ) 



where 



(12) 



E (x) = l-e 



is the function used to subtract continuum, s is the continuum threshold and 
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where M\ and M| are the Borel parameters in p and o channels. Since we assume that 
m a 700 MeV, which is close to the p meson mass, we will set M\ = Mf = 2M 2 and 
hence u = 1/2. 

It follows from Eq. (11) that in order to determine the coupling constant given in 
Eq.(12), one needs to know the parameters f p and X a . Leptonic decay constant of the p 
meson is known to have the value f p = 0.18 GeV from experimental result of the p 
decay width [19]. The residue X a is determined from sum rules obtained in [12] 
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where the last term describes the single instanton contribution, Ki are the modified Bessel 
functions and r is the instanton size for which we will use r = (0.6 GeV) -1 (see [20]). 

In further numerical analysis we use the following values for the input parameters: 
(m q qq) = (-0.82±0.1) x 10~ 4 GeV\ (a s G 2 ) = (0.038±0.011) GeV\ (a s (qq) 2 ) = (-0.18± 
0.1) x 10~ 3 GeV 4 (see for example [10]). For the mass of the scalar meson we use m a ~ 
700 MeV . The dependence of X a on Borel parameter M 2 at three values of the continuum 
threshold so — (1-4; 1.6; and 1.8) GeV 2 is presented in Fig. (1). Since the Borel parameter 
M 2 has no physical meaning, we must look for the so-called stability region where sum 
rules is practically independent of M 2 . We observe from Fig. (1) that the stability window 
M 2 lies in the interval 1.2 GeV 2 < M 2 < 1.4 GeV 2 and we get 



\ a = (0.2 ± 0.02) GeV 2 



(14) 



3 



where the error can be attributed to the variation of the continuum threshold sq, Borel 
parameter M 2 and the errors in the condensates (qq) and (G 2 ). It should be noted that if 
we had used the values of the input parameters given in [8] we get 

X a = (0.16 ±0.02) GeV 2 . 

This difference arises from the omitted instanton contribution and terms in perturbative 
contribution in [8]. 

Having the values of X a and f p , our next goal is calculating the coupling constant g pai 
from Eq. (12). As is obvious from Eq. (12) the main input parameters of the sum rule is 
photon wave functions. We shall make use of the following expressions for the photon wave 
functions [13, 15] 

(p(u) = 6u(l — u) , 
9i{u) = ~(l-u)(3-u) , 

and for the magnetic susceptibility we use x — —4.4 GeV~ 2 [21]. 

In Fig. (2) we present the dependence of the coupling constant g pal on the Borel parame- 
ter M 2 at three different values of the continuum threshold: so = 1-4 GeV 2 ; 1.6 GeV 2 ; 1.8 GeV 2 . 
It follows from this figure that for the choices s = 1A GeV 2 and s = 1.8 GeV 2 the varia- 
tion in the result is about ~ 10%, i.e., the coupling constant can be said to be practically 
independent of the continuum threshold Sq. Furthermore the coupling constant seems to 
be insensitive to the variation of the Borel parameter M 2 . Along these lines, we calculated 
the coupling constant g pal and our final result is value 

g pai = (2.2 ± 0.4) , (15) 

where all possible sources of uncertainties are taken into account, namely, errors coming 
from determination of X a , from variation of the continuum threshold s , Borel parameters 
M 2 , neglected twist -3 photon wave functions and errors in the values of the condensates. 

Finally, we would like to present a comparison of our prediction on the coupling constant 
g pul with the existing results in literature. Traditional 3-point QCD sum rules analysis 
predicts g pa ^ = (3.2 ±0.6) [8]. Our result is about 50% lower than this result, which can be 
attributed to the difference in the values of X a . More over our result is slightly lower than 
the result obtained from effective Lagrangian approach, which predicts g p(T1 = 2.71 [7]. 
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Figure captions 



Fig. (1) The dependence of \ a on Borel parameter M 2 at three different values of the 
continuum threshold s = 1.4 GeV 2 ; 1.6 GeV 2 and s = 1.8 GeV 2 . 

Fig. (2) The dependence of the g pal coupling constant on Borel parameter M 2 at three 
different values of the continuum threshold so = 1.4 GeV 2 ; 1.6 GeV 2 and sq = 1.8 GeV 2 . 
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Figure 1 
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Figure 2 



